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Diagrammatic analysis for normal state of Hubbard model proposed in our previous paper' 1 ' is 
generalized and used to investigate superconducting state of this model. We use the notion of 
charge quantum number to describe the irreducible Green's function of the superconducting state. 
As in the previous paper we introduce the notion of tunneling Green's function and of its mass 
operator. This last quantity turns out to be equal to correlation function of the system. We proved 
the existence of exact relation between renormalized one-particle propagator and thermodynamic 
potential which includes integration over auxiliary interaction constant. The notion of skeleton 
diagrams of propagator and vacuum kinds were introduced. These diagrams are constructed from 
irreducible Green's functions and tunneling lines. Identity of this functional to the thermodynamic 
potential has been proved and the stationarity with respect to variation of the mass operator has 
been demonstrated. 

PACS numbers: 71.27.+a, 71.10.Fd 



I. INTRODUCTION 

The present paper generalizes our previous work' 1 ! on 
diagrammatic analysis of the normal state of the Hub- 
bard model' 2-4 ! to the superconducting state. 

Now we shall assume the existence of pairing 
of charge carriers and non-zero Bogolyubov quasi- 
averages' 5 ' and, consequently, of the Gor'kov anomalous 
Green's functions' 6 '. 

The main property of the Hubbard model consists in 
the existence of strong electron correlations and, as a 
result, of the new diagrammatic elements with the struc- 
ture of Kubo cumulants and named by us as irreducible 
Green's functions. These functions describe the main 
charge, spin and pairing fluctuations of the system. 

The new diagram technique for such strongly cor- 
related systems has been developed in our earlier 
papers[7 — 17]. This diagram technique uses the algebra 
of Fermi operators and relies on the generalized Wick 
theorem which contains, apart from usual Feynman con- 
tributions, additional irreducible structures. These struc- 
tures are the main elements of the diagrams. 

In superconducting state, unlike the normal one, the ir- 
reducible Green's functions can contain any even number 
of fermion creation and annihilation operators, whereas 
in normal state the number of both kinds is equal. There- 
fore we need an automatic mathematical mechanism 
which takes into account all the possibilities to consider 
the interference of the particles and holes in the super- 
conducting state. 

With this purpose we use the notion of charge quantum 
number, introduced by us in [7] and called a— number, 
which has two values a = ±1 according to the definition 



Were C is a Fermion annihilation operator. In this new 
representation the tunneling part of the Hubbard Hamil- 
tonian can be rewritten in the form 

ff ' = EE *(*' - (2) 
- \ E E E ^'-^,0^, 

a — —1,1 CT rj) ' 

with the definition of the tunneling matrix elements 

tx{lt'-^) = t(-j^'--#) 
t- l (-tf--&) = t{£-lP) (3) 
t(t = 0) = 0. 

In this charge quantum number representation the oper- 
ator H' has an additional multiple a for every vertex of 
the diagrams and additional summation over a. All the 
Green's functions depend of this number. 

In interaction representation operator H' has a form 

H '( T ) = ^E E (4) 
aa -US' 

x G^(r + aO+)C£ a (T), 

The main part of the Hubbard Hamiltonian 

H = H° + H', 
H° = (5) 

i 

contains the local part H°, where fi is the chemical po- 
tential and U is the Coulomb repulsion of the electrons. 
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This interaction is considered as a main parameter of the 
model and is taken into account in zero approximation of 
our theory. The operator H' describes electron hopping 
between lattice sites of the crystal and is considered as a 
perturbation. 

We shall use the grand canonical partition function in 
our thermodynamic perturbation theory. 

The paper is organized in the following way. In section 
II we determine the definition of one-particle Matsubara 
Green's functions by using a representation and develop 
the diagrammatic theory in the strong coupling limit. 

In section III we establish relation between the full 
thermodynamic potential and the renormalized one- 
particle Green's function in the presence of additional 
integration over auxiliary constant of interaction A and 
prove the stationarity theorem both for a special func- 
tional consisting of skeleton diagrams and for a renormal- 
ized thermodynamic potential shown to be its equivalent. 

II. DIAGRAMMATIC THEORY 

We shall use the following definition of the Matsubara 
Green's functions in the interaction representation 

G aa '(x\x') = (tC^(t)CZ?;Jt')U(P)) C o , (6) 

where x stands for (af,CT, r), index c of (■■■}q means the 
connected part of the diagrams and (...} means thermal 

average with zero order partition function ' gH ,i . 

We use the series expansion for the evolution operator 
U((3) with some generalization because we introduce the 
auxiliary constant of interaction A and use XH' instead 
H': 



Uxifi) = Tcxp(-A J H'(r)dr), (7) 
o 

with T as the chronological operator. At the last stage 
of calculation this constant A will be put equal to 1. 

The correspondence between definition (6) and usual 
one I 12 ' is the following: 

G\\x\x>) = -(TC^ a {T)C^Jr')U x (fi))l 
= G^,(^,t\^',t'), 

G\-\x\x') = -(TC^(r)C^Jr')Um)l 

= F^,(^,r\^',r% (8) 
G^ixW) = -(TC^(r)C^Jr')Um) C 

G?'-\xW) = -GY{x'\x). 

As a result of application of the Generalized Wick The- 
orem we obtain for propagator (6) the diagrammatic 
contributions depicted on the Fig.l In superconducting 



state, unlike the normal state, the propagator lines do 
not contain arrows which determine the processes of cre- 
ation and annihilation of electrons because indices a can 
take two values a = ±1 and every vertex of the diagram 
describes different possibilities. 

In Fig.l the diagram a) is the zero order propagator, 
the diagram b) and more complicated diagrams of such 
kind are of chain type. They correspond to the contribu- 
tion of the ordinary Wick theorem and give the Hubbard 
I approximation. The contributions of the diagrams c) 
and d) of Fig.l are 

c) = ^TC|,(r)C t 7(r 1 + ai 0+)^(r 1 )C_^'(r')^ r 
x axt ai { ~t' - ~t), 

d) = ^TC^(r)Cy, 1 (n)C|(7 2 )C^/(r')^ , 
x ait ai (f' - ~f)a 2 t a2 - ~f) 

where (...)Q r means the irreducible two-particle Green's 
function! 2-5 ! and summation or integration is under- 
stood here and below when two repeated indices are 
present. Spin index has been omitted for simplicity. 
In the diagram c) the equality of lattice sites indices 

af = T > ' = T > =af'is assumed and in diagram d) 

~£ = ~t = ~i = it' . The diagrams Fig.l c), d) and e) 
contain irreducible two-particles Green's functions, de- 
picted as the rectangles. In higher orders of perturba- 
tion theory more complicated many-particle irreducible 

Green's functions Gr»°^ ir [l, 2, n] appear. These func- 
tions are local, i.e. with equal lattice site indices. There- 
fore the diagram c) in Fig.l can be dropped since it con- 
tains a vanishing matrix element, t(x — x) = 0. The pro- 
cess of renormalization of the tunneling amplitude shown 
in the diagrams c) and d) leads to the replacement of the 
bare tunneling matrix element at a (x — it) in c) by a 
renormalized quantity T a a (x'\x). This process is deter- 
mined by the equation 

T^{x'\x) = T^ a Vk)+^f Q VM 0) 
x G^( Xl \x 2 )T£l^(x 2 \x), 

where 

T^'Vlz) = 8 aa ,at a {lt' - lt)8{r' -r- aO+)<W, 

(10) 

and G" 1 * 2 is the full one-particle propagator. The quan- 
tity T a a is shown in the diagrams as a double dashed 
line. 

We then introduce the notion of correlation function 
A aa (x\x') which is the infinite sum of strongly connected 
parts of propagator's diagrams. If we now omit from 
these diagrams all those contained in the process of renor- 
malization of the tunneling matrix element, we obtain 
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FIG. 1: The examples of the first orders perturbation theory diagrams for propagator. Solid thin lines depict zero order 
one-particle Green's functions and rectangles depict two- and four-particle irreducible Green's functions. Thin dashed lines 
correspond to tunneling matrix elements. Double solid line corresponds to renormalized propagator. 



the skeleton diagrams for correlation function. In such 
skeleton diagrams we replace thin dashed lines by dou- 
ble dashed lines and obtain the definition of A aa (x\x') 
shown in the Fig. 2 

There are two kinds of A dependence in the diagrams 
of Fig. 2. One is conditioned by dependence of T" Q and 
the second is determined by A being an explicit pre- 
factor in the diagrams. In Hubbard I approximation 
only the free propagator line is taken into account. All 
the contributions of Fig. 2 except the last one are local 
and their Fourier representation is independent of mo- 
mentum. Only these diagrams are taken into account 
in Dynamical Mean Field Theory! 18 !. The last diagram 
of Fig. 2 has the Fourier representation which depends of 
momentum. 

As a result of diagrammatic analysis we can formu- 
late the Dyson-type equation for full one-particle Green's 
function (x = "a^,r): 

j3 

aia2 

x KZl{x\xi)Ti%^{xi\x 2 )G^{x 2 \x'). 

(11) 

This equation can be written in the operator form: 

G= {l-Kf a y 1 K = l{l-f a K)- 1 (12) 

Using equations (9) and (12) we obtain the Dyson 
equation for the tunneling Green's function 

f = T°(1- AT )- 1 = {l-T°l)- l T°, (13) 

where the correlation function A has the role of mass 
operator for the renormalized tunneling Green's function. 

In Appendix A we demonstrate the equivalence of the 
equation (11) to usual ^ representation of superconduct- 
ing Green's functions. 



III. THERMODYNAMIC POTENTIAL 
DIAGRAMS 

The thermodynamic potential of the system is deter- 
mined by the connected part of the mean value of evolu- 
tion operator: 

F(\)=F -±(U X ((3)) C , (14) 

with A equal to 1. 

In the Fig. 3 are depicted the first order diagrams for 

The notations in Fig. 3 are n = (jf,T n ) and n' — 
(lt',T n + a n + ). 

The first three diagrams in Fig. 3 are of chain type and 
correspond to the Hubbard I approximation. The next 
diagrams contain the rectangles which represent our ir- 
reducible Green's functions. Indeed, some of these dia- 
grams are equal to zero when the dashed lines are self- 
closed by virtue of the relation i(0) = 0. However, when 
these dashed lines are replaced by renormalized quanti- 
ties T\ their contributions are different from zero and 
should be retained. Such renormalized tunneling quanti- 
ties will be used in the next part of the paper. The contri- 
butions of the fifth and eighth diagrams on the right-hand 
side of Fig. 3 are 

x ait ai (~t - l)a 2 t a2 (2 -2)a 3 t a3 (3 -if) 
x G ( ^ a2a3 (t,T 2 \f',T 3 ), 

+ 1 (TC^(ri)C^(r2)C^(r 3 )C^(n)Y 

x ait ai ( 1 - ~t)a 2 t a2 (2 - 2) 

x a 3 t a3 {3 - 3)a 4 i a4 (4 - 4), 
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FIG. 2: The skeleton diagrams for correlation function A aa (x\x'). The rectangles depict the many-particles irreducible Green's 
function. The double dashed lines depict the full tunneling Green function T" a (x\x'). 
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FIG. 3: The first orders of perturbation theory 
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contributions. 



respectively. 



and as the result we have 



Comparison of the diagrams of Fig.l for the G^ aa ' 
(n-th order of perturbation theory for the one-particle 
propagator) to the contributions of Fig. 2 for (U™ +1 (/3)) Q 
((n + l)-th order for evolution operator) allows us to es- 
tablish the following simple relation (n > 1): 



U 



(n+l) 



(/?) 



x G ^ rCl 



«i0 + ), 



1 1 



Ai 



o 



A 



T I 



(15) 



aiO + ) 



Taking into account equations (12) and (13) we obtain 



1 f d\x , 



(UxWWo = -o / -^Tr(\xT Xl A Xl 



(16) 
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Then from (14) and (16) it follows that 



As a consequence of these equations we have 



A 

F(X) = Fo + ^J ^Tr(A!f Al A Al ) 
o 

A 



where 

S A - AA A (18) 

has the role of mass operator for tunneling Green's func- 
tion T\. For them Dyson equation exists: 

f = f + T ST. (19) 
Equation (17) can be rewritten in the form 



A 



dF(\) 
dX 



2(3 



Tr(T\Ti\). 



(20) 



The equations (15) and (17) establish the relation be- 
tween the thermodynamic potential and renormalized 
one-particle propagator G\ or tunneling Green's func- 
tion T\. Both these quantities depend on auxiliary pa- 
rameter A which is integrated over. As have been proved 
by Luttinger and Ward' 19,20 ', for normal state of weakly 
correlated systems, it is possible to obtain another ex- 
pression for the thermodynamic potential without such 
additional integration. 

In our previous paper M, for the normal state of Hub- 
bard model, we have obtained such an equation in the 
form of special functional. We now consider its general- 
ization to the case of superconductivity. For this purpose 
we introduce the functional 



Y(X)=Y 1 (X) + Y'(X), 



(21) 



where 



Ki(A) = --Tr{ln(AT°A A - 1) + T A AA A }, (22) 

and Y'(X) is the functional constructed from skeleton di- 
agrams depicted on Fig. 4. 

From Fig. 4 and Fig. 2 it is possible to obtain the rela- 
tion 



5T« a \x\x') 2' 



(23) 



Now we shall take into account the following functional 
derivatives based on the equation (12) and (13): 

f— — -Tr(ln(f °AA - 1)) = -Tr(f - ^ , ), 

ST« a '(x\x') V V " y 6T« a '(x\x'" 



5T% a '(x\x') 



—Tr(TAX) = A A" a (x'\x) 



+Tr(T x 



5A> 



5T« a (x\x 



7^ 



(24) 
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8T% a '{x\x') 



and 



— Tr{ln(AT°A A — 1) + AT A A A } = XA a x a (x'\x), 



(25) 



5T« a (x\x' 



(26) 



With the functional derivative of Y'(X) given in (23) we 
obtain the stationarity property of the functional l^A): 



SY(X) 



6T£ a \x\x') 







(27) 



Using the definition (18) of the mass operator S A we 
can rewrite the functional Y\ (A) in the form 

li(A) = -iTr{ln(f°£ A - 1) +f A S A }, (28) 

and prove the second form of stationarity property 
6Y(X) 



<5£> 



= 0. 



(29) 



To demonstrate this equation it is sufficient to use the 
Dyson equation (19) in the form 



T° = f - 1 + E A , 
and the derivatives: 



(30) 



S(T-yP'(y\y>) 



<5(E A ) Q «' {x\x') 



— — SapSa'p'SxySx/yl, 



6(£ x )™'(x\x>) 
5 



S(Z x )<*"'(x\x') 



(f x r'P'(x'\y')(f x )^(y\x), 



Tr(T x X x ) = (T x ) aa (x'\xX31) 



+ (T x E x T x ) aa (x / \x), 



5{Z x ) aa ' {x\x' 
Therefore we have 
SYiW 



Tr{ln(T°X x -l)} = -(T x ) aa (x'\x). 



= --(T x Y, x T x ) aa (x'\x), (32) 
5(Z x ) aa '(x\x') 2 v a a x } v 1 



and 



6Y'(X) 



5Y'(X) 5(T x )W(y\ y ') 



5(Z x )™'(x\x>) 5(T x )W'(y\y>) 5(Z X )<*°' (xtf) 

= hf x E x f x r' a (x'\x), (33) 



FIG. 4: The skeleton diagrams for functional Y'(X). The rectangles depict the irreducible Green's functions. The double 
dashed lines depict the tunneling renormalized Green's functions T aa (x\x ). 



where the usual convention about summation over the 
repeated indices has been adopted. 

As a result we obtain the stationarity property (29) of 
the functional Y(X) versus the change of the mass oper- 
ator This mass operator for A = 1 coincides with 
correlation function of our strongly correlated model. 

Now it is necessary to find a relation between the ther- 
modynamic potential F(X) and the functional Y(X). 

Consider first the value of the derivative ^jj~< The 
A dependence of the functional Y(X) is of two kinds: 
through Y,\ and also explicit through the factors A™ in 
front of the skeleton diagrams for the functional Y'(X). 

Due the stationarity property (29) we obtain 



dY(X) _ 5Y{X)dZ x dY(X). 



dX 



5Sa dX 
dY(X) 
~~dX~ 



dX 
dY'(X) 
dX 



(34) 



Here we took into account that the Yi(A) part of func- 
tional Y(X) (see equations (21) and (28)) does not ex- 
plicitly dependent on A. 

By using the definitions of Y' (A) (see Fig. 4) and of A\ 
(see Fig. 2) it is easy to establish the property: 



dY(X) dY'(X) 1 , ^ -> , 

From the equations (20) and (35) we have 
dY{X) 



X- 



dX 



—Tr(T\Y<\) = A — , 



and we therefore obtain 



F(X) = Y(X)+F , 



(35) 



(36) 



(37) 



since for A = the perturbation is absent Y(X = 0) = 
and F(X = 0) = Fq. Now we set A = 1 and obtain 



F = F +Y(1), 
with the stationarity property 
SF 



as 



= o. 



(38) 



(39) 



IV. CONCLUSIONS 

We have further developed the diagrammatic theory 
proposed for strongly correlated systems many years ago 
to establish the stationarity property of the thermody- 
namic potential in the superconducting state of the Hub- 
bard model. 

First, we have introduced the notion of charge quan- 
tum number which gives the possibility to consider the 
presence of irreducible Green's functions with an arbi- 
trary number of creation or annihilation Fermi- operators 
in superconducting state. 

We have introduced the notion of tunneling Green's 
function and its mass operator, which turns out to be 
equal to the correlation function of the fermion system. 

We have proven the existence of the Dyson equation for 
this function and establish the exact relation between the 
thermodynamic potential and renormalized one-particle 
propagator. This relation contains an additional integra- 
tion over the auxiliary constant of interaction A. 

We have constructed a special functional based on the 
skeleton diagrams for the propagator and for the evolu- 
tion operator which contain the irreducible Green's func- 
tions and full tunneling Green's functions. 

We have proven the existence of the stationarity prop- 
erty of this functional and establish its relation with ther- 
modynamic potential. 
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It is important to emphasize that there is a close simi- 
larity between our results obtained for two different mod- 
els of strongly correlated systems such as Periodic Ander- 
son Model (PAM) and the Hubbard Model (HM). From 
comparison of the results obtained for the PAM (see pa- 
per [17]) and the results of the present paper for the 
HM the topological coincidence of the diagrams for both 
models has been revealed. 

For example the skeleton diagrams of Fig. 3 of paper 
[17], obtained for A functional of PAM topologically co- 
incide with the skeleton diagrams of our Fig. 2 for the 
same functional, but of quite a different model. In or- 
der to obtain a complete coincidence, it is necessary to 
replace the full Green's function G c (iuj) of conduction 
electrons of PAM by the renormalized tunneling Green's 
function T{iuj) of the HM. 

The same similarity exists between other functionals 
of these models. For example, comparison of the skele- 
ton diagrams of Fig. 10 of paper [17] with the diagrams 
of Fig. 4 of the present paper reveals the full coinci- 
dence upon replacement of the Green's functions G c by 
T. This comparison allows us to conclude that from 
the thermodynamic point of view the Periodic Ander- 



son Model can be reduced to the Hubbard Model if we 
replace the Green's function of the conduction electrons 
of PAM subsystem by tunneling Green's function of hop- 
ping electrons of HM. This equivalence is irrelevant for 
the kinetic properties of PAM. 

We also note that the skeleton representation of our 
functional allows to select the local irreducible Green's 
functions as can be seen from Figs. 2 and 10. These 
quantities contain only fluctuations in time, unlike the 
non local ones which include both fluctuations in time 
and space. The coefficients of local diagrams (see Fig. 
10) vary with the order of perturbation theory as 2 „i lrt , 
for n > 1. 

Only such local diagrams are relevant for DMFT, so 
that one can attempt to carry out the summation of this 
class of diagrams. 
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Appendix A: Gor'kov-Nambu representation A = ±1 we obtain (k = (Jc, iu> n )): 

We consider equation (11) in Fourier representation. 
By inserting specific values of charge quantum number 



G£>) = A*<» + A&(fc)ei(t)G^(*0 - Ai:ZUk)e-i(t)GZl'lAk), (Al) 
Gl'-_Uk) = Al : zl(k) + A^(fc) ei (t)G^:\,(fc) + Ai;: 1 CTi (fc) e _ 1 (t)G^,_ (T (-fc), (A2) 

G:i;Ufc) = A:£,(fc) + A:^^ (as) 

- A^, - AI^ i (fc) £l (t)G^ (T ,(fc) + A 1 ji li _ <r (-fc)e_ 1 (t)G^, i _ tri (-fc). (A4) 

I 



Here 



e 1 (t) = e(t), c _ 1 (t) = 6(-t), £ (t) = l^^)e^, £e(t) = 0, G^ 1 ;" 1 ^) = -G^(-fc). (A5) 



r 



Assuming that the system is in a paramagnetic state, (k) = Y^ a {k), (A6) 

that superconductivity has a singlet character and using 
the definitions (8) together with the additional ones: 

A^(fc) = K a {k),k 1 a - 1 (k) = Y a ^(k), we obtain the following results: 
I 



r (1A _ A CT (fc)(l - e(-t)A w (-k)) - e(-t)Y aW (k)Y^(k) _ Y aW (k) v _ Y Wa (k) 

Ga{k) - iAk) ' ^^'-^M' ™ [) ~^w [ ' 

d a (k) = (1 - e(t)A CT (fc))(l - e(-t)A-(-fc)) + e(t)e(-t)Y^(k)Y^(k), 



which coincide with those found in the papers^ 9 ' 10 ' . By using equation (9) we can obtain 

In spinor representation the system of equations 
(A1-A7) has the form 



G = A + AeG, (A8) 



_ c(¥)(l - e(-¥)A ? (-fc)) 



were a , dJk) 

T£\k) = -e(t)e(-t)FMk), (A10) 

x _ / K{k) Y aW (k) \ ? _(e{k) V 



Y^ a (k) -A w (~k))^ V -c(-fc) 



